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ABSTRACT 

The Nuclear Compton Telescope (NCT) is a balloon-borne Compton telescope designed for the study 
of astrophysical sources in the soft gamma-ray regime (200 keV-20 MeV). NCT's ten high-purity 
germanium crossed-strip detectors measure the deposited energies and three-dimensional positions of 
gamma-ray interactions in the sensitive volume, and this information is used to restrict the initial 
photon to a circle on the sky using the Compton scatter technique. Thus NCT is able to perform 
spectroscopy, imaging, and polarization analysis on soft gamma-ray sources. NCT is one of the next 
generation of Compton telescopes — so-called compact Compton telescopes (CCTs) — which can 
achieve effective areas comparable to COMPTEL's with an instrument that is a fraction of the size. 
The Crab Nebula was the primary target for the second flight of the NCT instrument, which occurred 
on 17-18 May 2009 in Fort Sumner, New Mexico. Analysis of 29.3 ks of data from the flight reveals 
an image of the Crab at a significance of 4<j. This is the first reported detection of an astrophysical 
source by a CCT. 

Subject headings: gamma rays: general — ISM: individual (Crab Nebula) — Techniques: imaging 
spectroscopy — Methods: data analysis — Telescopes — Balloons 



1. INTRODUCTION 

Soft gamma rays (200 keV-20 MeV) are important 
probes of the most violent and extreme processes in the 
cosmos. Gamma rays are produced by non-thermal pro- 
cesses in such disparate objects as neutron stars, X-ray 
binaries, and Active Galactic Nuclei (AGN), and they 
also result from the decays of many radioactive nuclei 
that are released by supernova explosions, such as 26 Al, 
60 Fe, and 44 Ti. The penetrating nature of gamma rays 
allows the astrophysicist to probe deep within these often 
obscured systems and make unique and complementary 
observations of their gravitational fields, magnetic fields, 
and nuclear reactions. 

Over the last four decades, various types of telescopes 
have been developed to detect and image soft gamma- 
rays. One prominent technology is the Compton tele- 
scope, which exploits the Compton effect to perform di- 
rect imaging of gamma-ray photons. The Compton effect 
occurs when a gamma-ray photon incoherently scatters 
with an electron and is deflected by an angle (p. The 
Compton scatter formula relates the energy of the scat- 
tered gamma-ray photon E 7l and the incident photon 
energy E lQ to the scatter angle ip: 
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wher e m e c 2 = 511 keV is the rest mass energy of the elec- 
tron (|Comptonl [19 23) . The basic principle of Compton 
imaging is to measure the energies deposited by at least 
one Compton scatter and employ Equation Q] to recon- 
struct the direction of the incoming photon. If both the 
direction of the scattered photon and cos cp can be mea- 
sured, then the photon's direction can be constrained to 
a circle on the sky (the "Compton circle"). The Comp- 
ton effect is also intrinsically sensitive to the polariza- 
tion of the i ncident photon through the Klein-Nishina 
cross section (Kl ein fc Nishinalll929[ ). enabling a Comp- 
ton telescope to potentially measure the polarization of 
a gamma-ray source. 

The first satellite-borne Compton telescope was 
COMPTE L on board the Compt on Gamma Ray Ob- 
servatory (jSchonfelder et al.l Il993h , which consisted of 
two layers of scintillator detectors spaced by 1.5 me- 
ters. Photons were required to scatter in the top layer 
and be fully absorbed in the bottom layer. COMPTEL 
succeeded in making measurements and images of the 
gamma-ray regime to unprecedented sensitivity, includ- 
ing t he first detection o f 44 Ti from a supernova rem- 
nant (llvudin et al.lll994f) and a m ap of the galactic 26 Al 
distribution (jOberlack et~a"T1 Il996[ ) . 

To build on the successes of COMPTEL, the next 
generation Compton telescopes utilize a compact de- 
sign, which requires detectors with both fine position 
and fine energy resolution to track individual gamma- 
ray interactions within the detector volume. Com- 
pact Compton Telescopes (CCTs) can achieve effec- 
tive areas comparable to COMPTEL's with an in- 
strum ent that is a fraction of the size (Bo ggs fc Jeanl 
12001? ) . The first CCT was LXeGRIT, which uses a 
single liquid xenon time projection chamber (TPC) 
to measure interaction positions with high spatial res- 
olution (^3 mm) and energy r esolution («10% at 
1 MeV) (lAprile et al.lll993aL l2000h . Various CCTs are 
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terials (e.g., iKurfess et all (12007ft: 


O'Neill et all (12003): 


iBloser et al.| (|2002h: Takeda et al. 


(20071): lUeno et all 



IKurfess et all ([20031 )). 

In this paper we will briefly describe the NCT instru- 
ment (Section EJ), then describe the event selections and 
present the first image of the Crab (Sections [3] and 2]). 
We discuss two analysis tools: Monte Carlo simulations 
of our observation of the Crab and the Angular Resolu- 
tion Measure (ARM) for a Compton telescope (Sections [5] 
and |B|). Using these tools, we determine the significance 
of the Crab detection using the total counts, an analysis 
of the spectrum, and an analysis of the ARM histogram 
(Sections [THOj). Polarization analysis of the Crab was not 
performed due to the lack of sufficient counts. 

2. THE NCT INSTRUMENT 

The NCT instru ment (jBoggs et al.ll200lL I2003L 120041 : 
iChang et all I2007D is a balloon-borne CCT that uti- 
lizes high-purity crossed-strip germanium detectors for 
Compton imaging. To date, NCT has flown on two 
conventional balloon flights. The first was a proto- 
type flight in 2005 that succeeded in measuring the soft 
gamma-ray at mospheric background and galact i c an- 
ticenter regi on (jCoburn et al.l 120051 : iBowen et all 120071 : 



lBowenll2009l ). The full 10-GeD ve rsion of NCT was flow n 
on its second flight in May 2009 (jBandstra et al.ll2009f ). 
and the first results are presented here. 

The heart of NCT is its ten cr ossed-strip high- 
purity germanium detect ors (GeDs) (j Amman fc Lukd 
120001 : lAmman et all 12007ft . Each GeD has the ability 
to record in three dimensions the location of each in- 
dividual photon interaction. A single GeD measures 
8 cm x 8 cm x 1.5 cm and has 37 strips on each 
side, with a strip pitch of 2 mm. NCT is designed 
such that when a gamma ray undergoes at least one 
Compton scatter and is fully absorbed inside the de- 
tectors, the individual scatter positions and energy de- 
posits can be measu red. Compton Kinematic Discrim- 
inati on algorithms (f Aprile et al. 1993b; Bog gs fc Jeanl 
120001 : lOberlack et aLll2Q0Qt iKroeger et al.ll2QQ2h are used 
to find the most likely scattering order, and then the 
total gamma-ray energy can be reconstructed and the 
initial direction can be traced to a ring on the sky. Thus 
NCT works as both a spectrometer and an imager, and is 
intrinsically sensitive to polarization because it exploits 
Compton scattering. 

The entire set of detectors and their cryostat are en- 
closed inside a well of anticoincidence BGO shields to 
reduce the Earth albedo and atmospheric backgrounds. 
The resulting overall field of view is primarily limited 
by the BGO shields to ~3.2 sr. The entire instrument 
and readout electronics are mounted in a pointed, au- 
tonomous balloon gondola. Further deta i ls of th e gon- 
dola systems can be found in Bellm et al.l (|2009af ). 

Data from the detectors are proces sed on-board by con- 
ventional surface-mo unt electronics (jCoburn et al .1120041 : 
iBandstra et al.ll2QQ7h . which return the energy and tim- 
ing of each strip above threshold during an event. The 
event information is processed on-board and telemetered 
to the ground. More details on the digital system, includ- 
ing flow charts of the d igital logic involve d in acquiring 
events, can be found in iHung et al.l (|2009[ ). 



The software analysis tools for NCT are built us- 
ing the Med ium Energy Gamma-ray A stronomy library 
(MEGAlib) (jZoglauer et al.ll200l l2008h. MEGAlib pro- 
vides utilities for simulations (Zoglau er et al.l 12009). ge- 
ometry mode ling, event reconstruc tion, and imaging re- 
construction (Zo glauer et al.l 12010). in addition to tools 
for building a custom analysis pipeline. The NCT 
da ta analysis pip e line is explained in further detail 
bv IBandstra et al.l (|2009f ). 

3. DATA SELECTIONS 

NCT was launched from the Columbia Scientific Bal- 
loon Facility (CSBF) in Fort Sumner, New Mexico 
(34.5°N, 104.2°W) at 1330 UT 5/17/2009. The flight 
was terminated at 0400 UT 5/19/2009 near Kingman, 
Arizona (34.9°N, 113. 7°W). Altitude was maintained 
between 35 and 40 km for the entire time the detectors 
were operational. The total flight time was 38.5 hours, 
with nine of the ten detectors on for a total of 22 hours 
at float. Approximately 33 ks of data were taken with 
the Crab above 40° elevation. A discu ssion of the flight 
can be found in IBandstra et al.l (j2QQ9h . 

Data from the flight were selected only when the Crab 
was above 40° elevation both to minimize atmospheric 
absorption and to sample the background (as will be dis- 
cussed in Sections [8] and |9|) . During the chosen obser- 
vation time, the average livetime for the instrument was 
nearly constant with an average value of 92%. Taking 
into account the livetime and ^1 ks of dead time due to 
data acquisition system and flight computer resets, the 
effective Crab observation time was 29.3 ks. 

3.1. Event Reconstruction 

Event reconstruction for events with three or more in- 
teractions was done using the Co mpton Kinema t ic Dis - 
crimination (CKD) algorithm of Bo ggs fc Jeanl (j2000f ). 
For events with TV interactions, there are N — 2 Comp- 
ton scatter angles that can be measured (i.e., at all but 
the first and last interaction points), and this allows for 
a figure-of-merit with N — 2 degrees of freedom to be 
calculated for each possible event ordering. The best 
interaction ordering of the TV! possible orderings is se- 
lected. In principle, the more interactions there are, the 
more information there is about the event and the bet- 
ter it can be reconstructed (or rejected) using CKD. A 
limit of seven interactions was chosen because the num- 
ber of events with more than seven interactions becomes 
very small (fewer than 1%) and the time to reconstruct 
a single event scales like N factorial. 

Two- interaction events (i.e., a Compton scatter and 
photoabsorption) lack any kinematic degrees of freedom 
and cannot be reconstructed using CKD. Their interac- 
tion ordering is often ambiguous, and a single-scatter re- 
construction technique must be used. After trying differ- 
ent techniques on the NCT calibration data, we found the 
technique that properly reconstructed the most events 
was to choose the ordering that maximized the product 
of the Klein-Nishina cross section for the implied scat- 
ter angle with the probability of a photoabsorption at or 
before the second interaction. 

3.2. Background Reduction 

The gamma-ray background at balloon altitudes was 
measured by the NCT prototype during its 2005 balloon 
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flight (|Bowenir2QQ9[ ). The background primarily consists 
of three kinds of background events: real photons com- 
ing from outer space and from cosmic-ray interactions 
in the atmosphere; cosmogenic particles such as protons 
and neutrons; and instrumental events coming from ra- 
dioactivities in the instrument. 

Simulations are used to estimate how many instrumen- 
tal events survive the event reconstruction process. Such 
comparisons to simulations have been done, and the re- 
sult is that instrumental events are effectively eliminated 
by CKD, and the background is dominated by real pho- 
tons - predominantly atmospheric continuum photons, 
cosmic photons, and atmospheric 5 11 keV photo ns; and 
particle interactions are very few ([Bowenl |2QQ9[ ) . The 
2005 measurements yielded good agreement with simu- 
lations, and these simulation inputs were used again to 
estimate the background during the 2009 flight. 

The atmospheric background continu um h a s bee n 
extensively stu di ed an d measured (e.g., iLind (|1975[) : 
iLing fc Gruberl fiS)Tfo IGraser fc Schonfelderl Mf% : 
iSchon felder e t all (fl977[ )). The continuum has an eleva- 
tion angle dependence — the specific intensity should be 
approximately constant and then rise sharply at the hor- 
izontal. The intensity in the horizontal direction should 
be larger than that from the vertical by a factor of ^2-4 
over the energy range 300-1000 keV. This sharp rise in 
the horizontal direction is another reason for requiring 
that the Crab be above 40° in elevation. 

In addition, the NCT's BGO shielding was designed so 
that it would have a very wide field of view away from 
the Earth, reducing background photons that are com- 
ing from the atmosphere below and from the sides. An 
Earth horizon cut (Table [T]) is also applied to further 
reduce background that is coming from the atmosphere 
below. However, background photons from the atmo- 
sphere above can only be reduced by shielding, with a 
corresponding reduced field of view. 

An energy cut of 300-495 keV and 520-1500 keV 
excluded a few other sources of background. The cut 
avoids the very strong atmospheric 511 keV background 
line. Also, below ~300 keV, discrepancies were noticed 
between the observed background continuum and sim- 
ulations, with NCT o bserving few er counts by a factor 
of two in that region (Chiu 2010). The comparison of 
simulations and effective area calibration data at higher 
energies reveals a deficit of measured photons when a 
source is placed in the plane of the detectors, presum- 
ably because the simulations underestimate the shield- 
ing effect of the detector guard rings. This effect is more 
pronounced for lower energy photons and likely leads to 
the measured deficit. This low-energy region (<300 keV) 
will be excluded from the analysis. 

3.3. Interaction Distance 

It also became clear from simulations and measured 
background rates that the other data cuts used for the 
flight data should be wide in order to attain the best 
sensitivity. We chose one major constraint — to limit 
the distance between the first and second interactions 
in the detectors — in order to improve the angular res- 
olution of the instrument and thus the signal-to-noise. 
This choice is effective because the position resolution 
of the NCT detectors is the the dominant con tribution 
to the angular resolution (Ba ndstra et al.l l2006). and the 



Table 1 

Data cuts used on flight data 



Parameter 


Allowed range 


Crab elevation 


>40 u 


Photon energy 


300-495 keV, 520-1500 keV 


Number of interactions 


2-7 


Earth horizon 


Reject if >99% of Compton 
circle is below elevation of 0° 


Compton scatter angle (cp) 


no restriction (0°-180°) 


Minimum distance between 
first two interactions 


1.2 cm 


Minimum distance between 
any two interactions 


0.4 cm 



vector between the first two interactions sets the axis of 
the Compton circle. Therefore, requiring a longer initial 
scatter distance improves the angular resolution, but the 
cost is that the number of events is reduced, especially 
those with three or more interactions. Based on simu- 
lations of the Crab, setting a lower limit of 1.2 cm for 
the first interaction distance improves the angular width 
of the Crab (10° ARM FWHM to 7.4° ARM FWHM), 
while eliminating only 40% of the events within the cen- 
tral FWHM of the ARM distribution (see Section [6] for 
a definition and discussion of the ARM). The cut is also 
meant to improve the likelihood that an event is properly 
reconstructed by ensuring that the first two interactions 
are separated by multiple strips. 

The different cuts used in the analysis of the NCT flight 
data are summarized in Table HJ After applying these 
data selections, 65.8% of the remaining events had two 
interactions, 25.9% had three interactions, 6.6% had four 
interactions, and 1.7% had five or more. 

4. IMAGE OF THE CRAB NEBULA 

After applying the data cuts, an image was produced 
from the remaining data. Since NCT measures the 
location of the photon to a circle on the sky (not a 
point), Compton imaging techniques are used to pro- 
duce images. Compton imaging can be done through 
the straightforward backprojection of Compton circles 
onto the sky or with more sophisticated image process- 
ing. One imaging tool that has been used for NCT data is 
List-Mode Ma ximum Likelihood Expec tation Maximiza- 
tion (MLEM) (|Wilderman et al.l 119981). which h as been 
implemented in MEGAlib (jZodauer et al.ll2010f ). This 
method is iterative, and it uses a maximum likelihood 
statistic to refine the raw backprojection. The result is 
that point sources become clearer and other features are 
sharpened. 

We constructed an image from the NCT flight data 
using the cuts mentioned previously. We plotted a total 
of 289,128 event circles in a 80° x 80° backprojection, 
of which approximately 3,800 event circles (1.3%) are 
expected from the Crab based on simulations (see Sec- 
tion [5]). We performed a total of five iterations of the 
MLEM algorithm on the image, which is the approxi- 
mate number of iterations needed so that the FWHM 
of a point source image matches the 7.4° ARM FWHM. 
The result is that the Crab appears in the image (Fig. [I]), 
though it was not originally apparent in the backprojec- 
tion (Fig. [2]). The Sun was in the field of view but no 
source appears there, as expected. 
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The Crab image reveals a slight offset (~ 2°) in the 
position of the source. The significance of this offset 
can be evaluated in the following manner. Given the 
equivalent gaussian sigma of the ARM histogram for 
the Crab (7.4°/2.35 = 3.1°) and the signal-to-noise es- 
timate of 4cr (Section [7]), the one-sigma uncertainty of 
the i mage position c an be estimated as « 3.1°/4 = 0.8° 
(e.g., Bobrofj (|1986[ )). This statistical uncertainty, com- 
bined with the « 1° systematic uncertainty in the gon- 
dola aspect, means the « 2° offset in the source position 
is not statistically significant. 

5. SIMULATION OF THE CRAB OBSERVATION 

A simulation of the Crab Nebula spectrum was per- 
formed as a point of comparison for the NCT flight 
data. A detailed mass model for NCT has been devel- 
oped, along with code that applies the d etector response 
to the ideal simulate d detector hits ([Bandstra et al.l 
l2QQj IChiu et alJl2QQ9h . The MEGAlib simulation tool 
Cosima was used, which is based on the gamma-ray tele- 
scope libra ry MEGAlib and the M onte Carlo framework 
GEANT4 (jAgostinelli et all 12003( 1. A point source with 
a spectrum identical to the C rab Nebula spectrum mea- 
sured by INTEGRAL/SPI (jJourdain fc Roquesl 120091 ) 
was simulated at infinity along the z-axis of the detec- 
tors (i.e., the center of the field of view) for an equiva- 
lent real time of 100,000 seconds. The total Crab flux 
between 200 keV and 10 MeV was calculated to be 
2.194 x 10 -2 ph cm -2 s _1 . The total number of events be- 
tween 200 keV and 10 MeV that interacted in the detec- 
tors were 219,987, with 18,622 events passing the event 
cuts (8.5%). 

Compensation for atmospheric absorption was done af- 
ter the simulation by calculating the transmission prob- 
ability as a function of energy and multiplying the prob- 
ability with the spectrum of events that passed the data 
cuts. The more accurate way of accounting for atmo- 
spheric absorption is to perform Monte Carlo simulations 
with a layer of air of the appropriate column density, 
but a brief exploration of such simulations revealed that 
there was little difference. An average column density 
of 3.2 g cm -2 was calculated for the chosen observation 
time, and the Crab spectrum at 45°, 55°, 65°, and 75° 
was calculated (the Crab elevation peaks at 80° at tran- 
sit). An average spectrum was calculated, and it con- 
tained 71% of the total counts of the original spectrum. 
This transmission percentage was used as an approxima- 
tion for atmospheric absorption in the simulations. 

6. THE ANGULAR RESOLUTION MEASURE (ARM) 

An important tool that will be used in the data anal- 
ysis is the Compton telescope Angular Resolution Mea- 
sure (ARM). The ARM is a statistic that is used to 
meas ure the angular resolutio n of a Compton telescope 
(e.g JSchonfelder et al.l (|1 9931 ) ). An ARM histogram can 
be constructed when the source's location is known. The 
ARM histogram is a histogram of the following quantity: 

— V^geo Brecon (2) 

where cp geo is the scatter angle derived from the known 
source location and (/? rec0 n is the Compton scatter angle 
obtained from event reconstruction (i.e., the Compton 
scatter formula). For a single event, A(/?arm measures 
the deviation of the event's Compton circle from the true 



position of the source. A graphical explanation of this 
definition is shown in Fig. 03 The angular resolution of a 
Compton telescope is specified by the full-width at half 
maximum (FWHM) of the ARM histogram. 

The FWHM of the ARM histogram is affected by 
detector energy resolution, detector position resolution, 
and Doppler broadening. Doppler broadening is the de- 
viation from the ideal Compton scatter formula caused 
by scattering off of bound atomic electrons (not free 
electrons). It is a fundamental limit to the angu- 
lar resolution of a Compton telescope. This broaden- 
ing h as been extensively s t udied for different materi- 
als in iZodauer fc Kanbachl (|2003f ). Because of NCT's 
excellent energy resolution, the contribution from the 
energy resolution is not important. The position res- 
olution of the NCT detectors (especially the 2 mm 
strip p itch) is the the dom inant contribution to the 
ARM ( Bandstra et al.l I2006D . This is one reason the 
event cut on the first interaction distance greatly im- 
proves the ARM FWHM. 

The ARM is dominated by the distance between the 
first two interactions, because that vector forms the axis 
of the Compton circle. Subsequent scatters do not af- 
fect the ARM greatly; if anything, events with multiple 
Compton scatters have a narrower ARM. This can be 
understood by noting that in order to undergo further 
scatters, the photon must lose less energy in the first 
scatter. Thus the Compton scatter angle for these events 
must be on average smaller, reducing the uncertainty in 
the ARM. 

Since the position resolution of all the detectors is the 
same (2 mm GeD strip pitch, and 0.5 mm depth reso- 
lution), and the GeD energy resolutions are similar ( 2- 
4 keV FWHM at 662 keV), the ARM is uniform through- 
out the detectors. The ARM is predominantly affected 
by the geometry of the event. In fact, one can estimate 
NCT's ARM FWHM by using the requirement of at least 
1.2 cm between the first two scatters and the germanium 
detector strip pitch of 2 mm. The FWHM should be 
=9.5°, which is comparable to the 12.8° measured 
for the 662 keV calibration source and the 7.4° mea- 
sured for the Crab itself. By comparison, the contribu- 
tion to the ARM due to the energy resolution (assuming 
^3 keV FWHM) should be =0.3°. 

The ARM will be used to select events from the loca- 
tion of the Crab and test points. The ARM distribution 
can usually be approximated by a Lorentzian distribu- 
tion with FWHM T, and the sensitivity of data cuts is 
maximized by selecting events that are within 0.7r of the 
center. For the data cuts we are using, T = 7.4°, so we 
use an ARM cut of ±5.2° to maximize sensitivity. 

7. SIGNIFICANCE OF THE CRAB DETECTION 

To measure the significance of the Crab, events that 
were consistent with the Crab were extracted from the 
data. Events whose Compton circles passed within the 
optimal ARM cut of 5.2° around the position of the Crab 
were selected. 

An estimate of the background counts is needed to 
determine the significance. Because the Crab moves 
through NCT's wide field of view, we can select nearby 
points to estimate the background at the Crab's position. 
An added complication is that because photons are only 
localized to a circle, a photon from an off-source point 
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Figure 1. Image of the NCT flight data in the region of the Crab, with the positions of the Crab and the Sun marked. Five iterations 
of the MLEM algorithm have been performed. There is a clear enhancement in the vicinity of the Crab. The color scale is an arbitrary 
intensity scale. 



may be consistent with the Crab as well. For this analy- 
sis, we selected eight points that were on a circle 8 = 20° 
away from the Crab at evenly-spaced intervals. Events 
were selected from within 5.2° of these points. The an- 
gular separation of 6 = 20° was selected so that we could 
sample points that were close enough to the Crab that 
they have similar exposure times, and far enough away 
so that fewer Compton circles would be consistent with 
both the Crab and these points. There is a trade-off be- 
tween choosing an offset angle 5 so that the background 
estimate is reliable (i.e., 5 should be small) and so that 
the off-axis source contribution is small (i.e., S should 
be large). The value of 20° was chosen in this work as 
the largest angle such that the interval (—5, +5) could 
be sampled for the later ARM analysis (due to the 40° 
minimum elevation of the Crab). 

Table [2] shows the total number of events that were ex- 
tracted from the on-source and the off-source positions. 
The results of the same analysis on the Crab simulations 
(i.e., no background) are also shown in the table. The 
number of excess counts observed at the position of the 
Crab agrees with the simulations to within statistical er- 
ror. 

The significance of the excess source counts can be es- 



Table 2 

The number of events with Compton circles within 5.2° of the 
Crab 



Position 


Counts (Flight 
Data) 


Counts (Scaled 
Simulation of Crab) 


(*o,6o) 


29,808 


1101 




29,041 


481 


Eight 
sample 
points 20° 
away 


28,708 
29,453 
29,462 
29,124 
28,917 
29,093 


459 
444 
453 
482 
457 
432 




29,329 


447 


Average 
off-source 


29,141 


457 


Excess 
on-source 


667 


644 



Note. — The point (£ 0l b ) = (184.56° 
of the Crab in galactic coordinates. 



-5.78°) is the position 



timated in the following way. The naive approach is to 
assume that all of the measurements (off-source and on- 
source) are independent Poisson processes. Let the Crab 
point be labeled as C, and the off-source points labeled 
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Figure 2. Raw Compton circle backprojection of NCT flight data in galactic coordinates, with the position of the Crab and the Sun 
marked. A total of 289,128 event circles are plotted, and « 1.3% of that number are expected to come from the Crab. Notice that there is 
no clear source at the location of the Crab before image processing. The color scale is an arbitrary intensity scale. 

statistical uncertainty on the excess counts: 




1 

^excess = N C ~ iV off _ source = N C ~ ~ ^ JV< (3) 

i—1 



^excess Nc 



, = 182.9 



64 



2=1 



(5) 



This approach is not exact because the Compton event 
circles are likely to intersect multiple points, meaning Nc 
and Ni are not independent Poisson measurements. We 
identify the following independent Poisson processes that 
can be used to calculate a valid statistical uncertainty: 



Figure 3. Schematic of the definition of the Angular Resolution 
Measure (ARM). Notice that the ARM can be positive or negative: 
A(fi < 0, since its reconstructed Compton circle circumscribes the 
source position, while A(f2 > because its Compton circle does 
not circumscribe the source position. 



as i = 1 ... 8. Then this approach yields an approximate 



C + k : An event intersects k off-source points 

and not the Crab point 
C + k : An event intersects k off-source points 

and the Crab point 



(6) 



where k = ... 8. Let Ng +k and Nc+k be the counts 
measured for each Poisson process. These measurements 
are listed in Table [3] Then the number of excess Crab 
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Table 3 

Data used in the calculation of the signal-to-noise ratio of the 
Crab data 
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Total 


103,357 


29,808 




Flight: 


Simulation: 


Excess on-source 


667 


644 


Background noise 


163.4 


163.4 (using 


to 




flight value) 


Significance 


4.1cr 


3.9<r 
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Figure 4. Simulated spectrum for events within 5.2° of the Crab 
using data cuts given in the text. The averaged spectrum measured 
at points 20° away is also plotted. The difference of these two 
distributions is the signature we expect in the flight data after 
estimating the spectrum at nearby points. 



counts can be written in the following way: 

-^excess — C -^off— source 



(7) 



1 8 

2 = 1 

8 / 8 8 \ 

: Nc+k - g ( E kN c+* + E kN c+k I (9) 

fc=0 \k=l k=l ) 

8 8 

: Nc+o + E ( 1 " J) N c+k -\± kNcilQ) 



k=l v 7 k=l 

The factors of k in the off-source sum are needed because 
each event tallied by Nc+k or N^ +k is counted at k dif- 
ferent off-source points. The statistical uncertainty for 
the excess counts becomes 

Excess = N C+ + E f 1 " I) 2 N C + k + E (l) 2 N &ti 
k=l ^ ' k=l ^ ' 



163.4 



(12) 



The result of this calculation is that the Crab was de- 
tected at the 4.1cr level, which is within statistical error 
of the expected 3.9a from simulations. 

8. ANALYSIS OF THE CRAB SPECTRUM 

To examine the spectrum of the Crab, a spectrum is 
formed for the on-source data, and an average off-source 
spectrum is also calculated. Error bars are calculated 
for each spectral bin in the same manner described in 
Section [7] and Table 03 The results of this procedure for 
the simulated Crab alone are shown in Fig. HI 

The chi-squared statistic is calculated for the mea- 
surement assuming there is no source present (xo) ano ^ 
the measurement assuming the Crab spectrum as de- 
rived from simulations (Xcrab)- When combined with 
the number of degrees of freedom of the data (V), each 
chi-squared statistic yields a P- value, or the probability 
that random fluctuations alone would yield a chi-squared 
of equal or greater value. The P- value for the null hy- 
pothesis quantifies whether a source is present (Po <^ 1) 
or not (Pq ~ 1), and the P- value for the Crab quantifies 
whether the measured spectrum is statistically consistent 
with the simulation (Pcrab ~ 1) or not (Pcrab < !)• An 



400 

300 

^200 
o 

° 100 


-100 



— Crab (Sim) 
Crab (Flight) 



Xl =18.0 (u=6) 
%Lab =2-5(^=6) 



400 



600 



800 1000 
Energy (keV) 



1200 



1400 



Figure 5. Measured spectrum for the Crab plotted with the ex- 
pected excess from simulations. The spectrum is measured by sub- 
tracting the average spectrum from points 20° away from the Crab. 
Error bars are calculated using Poisson statistics. The measured 
spectrum is inconsistent with zero (2.7<r) but is consistent with the 
simulated Crab spectrum (Pcrab — 0.87). The bins are 200 keV 
wide. 

equivalent gaussian sigma n can be derived from a P- 
value by equating the chi-squared P- value with the com- 
plementary cumulative distribution function of a gaus- 
sian random variable: n = \/2 erf -1 (1 — P). 

Performing this analysis on the flight data (Fig. [5]) re- 
veals a spectrum that is inconsistent with a null result 
with a P- value of Pq = 0.0062 and no = 2.7<r, and the 
spectrum is consistent (Pcrab = 0.87 and ncrab = 0.2<j) 
with the simulations. To test the spectral analysis 
method, three other test points were chosen that had 
similar exposure during the flight: (160°, -10°), (170°, 
+10°) and (175°, +20°). Each of these points yields spec- 
trum consistent with the null hypothesis (the respective 
no values are 0.6<r, 1.0a, and 0.2<r), and mostly incon- 
sistent with the Crab spectrum (no = 2.3a, 1.4a, and 
3.3a). The results are shown in Fig. [6] 

9. ANALYSIS OF THE CRAB ARM 

The measured ARM histogram of the Crab should also 
match the simulated ARM histogram of the Crab in or- 
der to further confirm the detection. In order to see if 
there are excess counts in the ARM histogram at the 
location of the Crab, one must first estimate what the 
ARM should be in the absence of a source. To do this, 
we used the same eight 20° off-axis points that were used 
in the previous analyses. As mentioned earlier, the value 
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Figure 6. Measured spectrum at three different test points, with 
the expected excess for the actual Crab position shown. The spec- 
trum is measured by subtracting the average spectrum from eight 
points 20° away from each test point. The spectrum at the three 
test points is consistent with zero. 

of S = 20° was chosen in this work as the largest angle 
such that the full interval (—5, +5) could be sampled in 
this analysis (due to the 40° minimum elevation of the 
Crab). The ARM histogram was calculated for each of 
those points, and the average was constructed. 

An ARM histogram that is calculated at a point sep- 
arated from a point source by an angle S (20° in this 
case) will have an ARM histogram that no longer re- 
sembles a Lorentzian, but it is a mostly flat distribution 
between ±5 that decays away beyond that region. So the 
points we used to estimate the ARM have photons from 
the Crab in them, and not an insignificant number. We 
must account for the stray counts from the Crab in some 
way or risk underestimating the flux. An example of this 
off-source ARM contribution is shown in Fig. [7| where 
the on-source Crab simulation is shown with the 20° off- 
source Crab ARM. The difference between the on-source 
and off-source ARM histograms is the signature that we 
expect to see in the flight data when subtracting back- 
ground estimates that are an angle 5 away from the Crab. 

The difference between the ARM histogram at the 
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Figure 7. Simulated ARM distribution for the Crab as seen by 
NCT using the data cuts given in the text. The FWHM of the 
simulated ARM distribution is 7.4° ARM FWHM, as mentioned 
previously. The averaged ARM distribution at points 20° away 
is also plotted. The difference of these two distributions is the 
signature we expect in the flight data after estimating the ARM at 
nearby points. 
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Figure 8. Measured ARM excess for the Crab plotted, with the 
expected excess from Crab simulations shown. The spectrum is 
measured by subtracting the average spectrum from eight points 
20° away from the Crab. Error bars are calculated using Poisson 
statistics. The measured ARM excess is significantly inconsistent 
with zero counts (Prj = 0.00025, or 3.7a) and is compatible with 
the Crab simulation (Pcrab — 0.90). The signal-to-noise ratio of 
the central bin is 4.2<r. These results confirm the 4<r significance 
found in Section 

Crab position and the average ARM histogram of the 
off-source points was formed, and the error bars of each 
bin were calculated in the same manner described in Sec- 
tion [3 and Table 03 The result is shown in Fig. The 
measured ARM excess is significant (Prj = 0.00025, or 
3.7a). The ARM distribution is also very close to the 
simulated expectation (Pcrab = 0.90). In addition, the 
central bin, which has the largest number of counts, has a 
signal-to-noise ratio of 4.2a, similar to the overall source 
significance derived in Section [71 The expected simula- 
tion value for the central bin is 3.9a, so the simulations 
are in close agreement with the flight data. 

To test the ARM analysis method, three other test 
points were chosen that had similar exposure during the 
flight: (160°, -10°), (170°, +10°) and (175°, +20°). Each 
of these points is consistent with the null hypothesis (the 
maximum value of no is 1.6<j) and inconsistent with the 
Crab simulation (the minimum value of ncrab is 2.2a). 
The results are shown in Fig. [9] 

10. CONCLUSIONS 

The Crab Nebula was detected by NCT at a signifi- 
cance of at least 4a during 29.3 ks of observation time in 
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to date (jVedrenne et alJl2003l ). 

Though NCT should be sensitive to the polarization of 
gamma-ray sources, a polarization analysis of the Crab 
was not performed here due to the lack of sufficient 
statistics. Recent gamma-ray observations have mea- 
sured a high degree of polarization for the Crab Nebula: 
- 45% in the 200-800 keV (IBIS) a nd 100-1000 keV 
(SPI) regimes (jForotl 120081 : iDeanl 12008). Given this level 
of polarization, and the estim ated instrumental m odu- 
lation factor for NCT of 45% (jBellm et al.ll2009bD . the 
expected polarization signal from the Crab should be 
~ 20% of the total signal. Therefore, with no other 
changes to the NCT instrument, an increase in observa- 
tion time by a factor of ~ 25 would be needed to detect 
the presence or absence of polarization at 4cr significance, 
and even more to measure the degree of polarization it- 
self. Further collimation of the instrument could reduce 
this time, at the cost of a reduced field of view. 

The NCT instrument suffered extensive damage during 
a launch attempt in Australia on a subsequent balloon 
campaign. Once the instrument has been rebuilt, the 
focus of future balloon flights is to perform polarization 
analysis on point sources and to image extended sources 
such as the galactic 26 Al emission. 
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